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Abstract Invasion waves are a fundamental build- 
ing block of theoretical ecology. In this study we 
aim to take the first steps to link propagation fail- 
ure and fast acceleration of traveling waves to crit- 
ical transitions (or tipping points). The approach 
is based upon a detailed numerical study of var- 
ious versions of the Fisher-Kolmogorov-Petrovskii- 
Piscounov (FKPP) equation. The main motivation 
of this work is to contribute to the following ques- 
tion: how much information do statistics, collected 
by a stationary observer, contain about the speed 
and bifurcations of traveling waves? We suggest warn- 
ing signs based upon closeness to carrying capac- 
ity, second-order moments and transients of localized 
initial invasions. 

Keywords: Critical transitions, invasion waves, 
propagation failure, Fisher-KPP, FKPP, SPDE. 



1 Introduction 

The propagation of waves has been a central topic in 
spatial ecology for a long time. A primary motivation 
arises from fronts where a new species is introduced 
into an environment or an existing species consid- 
erably extends its habit. A classical example is the 
spread of muskrats in central Europe [77] . Other doc- 
umented examples are butterflies and bush crickets 
in the UK [75] and the cane toad invasion in Aus- 
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tralia [63 . Also bacterial growth [54"Il2"o] shows very 
similar spreading and wave phenomena. The refer- 
ences in [31,76 contain even more examples. 

From a theoretical perspective a first ground- 
breaking result is the modelling of invasion waves via 
reaction-diffusion equations by |22] and Kolmogorov, 
Petrovskii, Piscounov [33] (FKPP) who studied the 
partial differential equation (PDE) 



du 



d 2 u 

dx 2 



u(l 



(1) 



There are many different aspects that could be in- 
cluded in a reaction-diffusion model which are very 
interesting to match theory and experiment; see [3Tj 
I3"3"II2THI53"] . Nevertheless, the basic guiding principles 
obtained from simple models are still highly rele- 
vant. Here we shall restrict ourselves to the study of 
the following stochastic partial differential equation 
(SPDE) 



du 

at 



d 2 u 
dx 2 



f(u) + 'noise'. 



(2) 



For now, the reader may just think of the classical 
FKPP nonlinearity f(u) = u(l — u) and some noise 
process that vanishes at zero-population level; for 
more technical details see Sections |3]H] The detailed 
choices are discussed later. 

The main theme of this paper is the interplay 
between invasion waves and so-called critical tran- 
sitions [7TH5] . Basically, critical transitions (or tip- 
ping points) are drastic sudden changes in dynamical 
systems; for some background and details see Section 
[2] The first major question is whether {l])-© can un- 
dergo a 'critical transition'. We discuss this question 
from a more technical perspective in Section [2] On a 
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heuristic level, one may just consider a parameter in 
([2]) that is slowly varying. Suppose that there exists 
a wave with positive speed for some parameter range 
while the wave is stationary (or reverses direction) 
for another parameter range. Whether an invasion 
reaches a new habitat or not can have drastically 
different consequences so one probably would like to 
refer to this situation as a critical transition. 

Another case we shall consider in this paper is 
the situation where the wave speed becomes infinite 
at a special parameter value. Hence, a small param- 
eter variation can cause a dramatically accelerating 
invasion wave. 

The next step is to check whether early-warning 
signs for a critical transition exist. In this context, 
changes in vegetation patterns have been the main 
motivation recently [4"Tll34] . There are only a few 
studies on early-warning signs for spatial systems 
[141113) . In fact, early-warning signs for noisy waves 
generated by SPDEs have not been considered yet. 
This paper makes a first step in this direction. We 
focus on transitions for the wave speed (e.g. prop- 
agation failure) as it controls when and where an 
invasion front appears. Although some detailed mea- 
surements of waves are available [35 , 51 it is very dif- 
ficult to obtain precise global empirical information 
[3"T1 p. 92] about a wave. Here we restrict ourselves 
to a single spatial observation location i.e. records 
by a single stationary 'ecological observer' over a 
fixed time interval. The general idea that one may 
obtain spatial conclusions from local observations is 
not new |21j . However, our detailed comparative nu- 
merical study of several different variants of ([2]) with 
a focus on local early-warning signs seems to be a 
completely new direction; for more details on the nu- 
merical methods see Section[5] The main themes and 
results from the numerical studies are the following: 

(a) A description of the statistics for early-warning 
signs in SPDEs with wave propagation failure 
based upon closeness to carrying capacity and 
second-order moments. 

(b) A comparative study of (a) for different noise 
types (white, space-time white) and different mul- 
tiplicative noise nonlinearities (parametric, finite- 
system size, etc.). 

(c) Investigation of statistics near continuous wave 
speed transitions (and their unpredictability) for 
Allee effect nonlinearities in the deterministic part 
of the FKPP SPDE. 



(d) Suggestion of transient minima to analyze wave 
propagation failure and wave speed blow-up. 

Beyond the technical contributions we also try to 
link different methodologies. We combine approaches 
from biological invasions, critical transitions, Fisher- 
KPP (and Nagumo) waves, SPDEs and numerical 
methods. This approach should also be helpful to 
link several, mostly distinct, communities such as 
theoretical ecology, waves in theoretical physics and 
mathematical methods for SPDEs. 

The paper is organized as follows. In Sections [2][5] 
we give brief reviews of the essential facts required 
for the remaining part of the paper. Due to the in- 
terdisciplinary aspects, the brief reviews seem neces- 
sary. Readers familiar with all the background may 
forward to Section [5] where the multiplicative noise 
case for the FKPP SPDE and statistical warning 
signs are studied. The nonlinear noise case is consid- 
ered in Section [7] and the Allee effect in Section [8] 
Section [9] on transient phenomena and the influence 
of initial conditions concludes the main part of the 
paper. In Section[10]a number of generalizations and 
open problems are listed. 

2 Background - Critical Transitions 

A primary motivation to study critical transitions 
(or tipping points) arose from ecology, e.g. due to 
the theoretical work of Scheffer and co-workers [73l 
1721174] . Then it became clear from many distinct ap- 
plied problems 7l[ as well as from abstract mathe- 
matical considerations 15 1(1] that many features for 
early-warning signs are generic across many dynami- 
cal systems. Recent studies of laboratory [16,82, and 
full ecosystem [3] experiments re-inforced this view- 
point. 

Here we recall a few aspects of critical transitions 
for finite-dimensional systems relevant for this pa- 
per. Consider the pitchfork bifurcation normal form 
[HI p.282] 

— = w' = fiw + w 3 . for ujiEl,|ieI. (3) 
dt 

The homogeneous trivial branch {w — 0} consists of 
stable equilibria for /i < and unstable equilibria 
for fj, > since the linearized system around w = 
is W' = fiW with solution W(t) = !F(0)e^. The 
bifurcation at fj, = is sub-critical with two unstable 
branches {w — ±\/— //} for /1 < 0. 

Consider a slow parameter variation fj/ = e with 
< e « 1 and /i(0) < 0. Orbits near the homoge- 
neous branch will reach a neighborhood of (w, (i) — 
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(0, 0) and then jump away quickly indicating a criti- 
cal transition [45J Fig. 3(c)]. Before the jump the sys- 
tem is slow to recover from perturbations ('slowing 
down') for fi < since W(0)e^ -> W(0) as fi -> 
for fixed t. For a deterministic system, it is impossi- 
ble to measure the slowing-down effect once it starts 
tracking the homogeneous branch {w = 0} i.e. it is 
exponentially close to w = 0. However, for a stochas- 
tic version of ([3]) given by 



w = nil) + w + 'noise' 



(4) 



the random perturbations can constantly kick the 
system away from the trivial branch. Extracting statis- 
tics from these perturbations can make the slowing 
down effect measurable [7T1I45] . This is one motiva- 
tion to study stochastic traveling waves ©. 

An important question is which bifurcation points 
or quantitative transitions we would like to classify 
as critical transitions. In multiple time scale systems, 
such as (J4j) augmented with fi' = e, the classification 
of local bifurcation points is relatively straightfor- 
ward (46, Sec. 2-3]. The mathematical classification 
and early-warning signs from |45 p 46 | can be applied 
to many pattern-forming bifurcations in spatially ex- 
tended systems on bounded domains. One first de- 
rives the amplitude equations on the domain locally 
[12] , Only a discrete set of eigenvalues occurs [36l 
p. 2 10] and the usual local bifurcations for a finite 
number of eigenvalues passing through the imagi- 
nary axis can often be applied. 

For patterns on unbounded domains the situa- 
tion is less clear. We do not offer any solution to 
this problem and consider an example to illustrate 
the difficulties. Consider a traveling wave solution 
u(x, t) = u(x - ct), e.g. for (JTJ), with (x, i)elxR+ 
which satisfies u(x, 0) = 1 for x < and u(x, 0) = 
for x > 0. Imagine a habitat [sijiCa] C K with 
x\p > and define the mapping 



l(u,T) 



1 



\x 2 - Xi\ 



\u(x, T)\dx. 



If the invasion wave spreads towards x — oo (s > 0) 
and saturates at the carrying capacity u = 1 then 
there exists a finite time Tj such that for all T > Ti 
we have I(u,T) = 1. If a slow parameter variation 
causes the wave to become stationary (s = 0) then 
l{u,T) = for all T > 0. Although this indicates 
how one may define one possible critical transition 
scenario for waves, the situation is actually unclear 
since for fixed T > one may have T(u, T) = for 



s > and s = 0. This illustrates again that global 
definitions are intricate [46j Sec. 8]. 

For this paper we simply rely on the intuitive 
notion that the transition to a standing wave and 
also the transition to wave speed blow-up are impor- 
tant in the context of critical transitions and early- 
warning signs. 



3 Background - FKPP Equation(s) 

A more general version of the PDE ([1]) studied by 
Fisher 22 as well as by Kolmogorov, Petrovskii and 
Piscounov [43] is given by 



(5) 



for u — u(x,t), (x,t) e K x [0, oo). The parameter 
D > controls the diffusion and if f(u; fi) — fif(u) 
then fi > can be interpreted as a growth rate. 
Of course, an initial condition has to be specified. 
Often one considers u{x, t — 0) either with compact 
support localized near x = or an initial condition 
with Gaussian decay. The localized initial condition 
for a population u to appear in a new environment 
is not only a mathematical simplification but does 
occur under realistic conditions e.g. due to global 
long-range transportation networks [44) . 

The nonlinearity / : R 2 — > K represents growth 
and saturation effects and is required to satisfy the 
conditions 

/(0 ;At ) = 0, /'(0;/i)>0, 
/(l;/i) = ! /'(l;/i)<0. 

The classical example is logistic growth /(«; fi) = 
fjbu{l — u). In this case one may rescale t n- t/fi, 
x i y x\J D/fj, to obtain from ([5]) the FKPP equation 



du d 2 u 



(6) 



Initially, the FKPP equation © modeled the spread 
of genes in a population but it has since become a 
paradigmatic model for populations dispersing un- 
der the influence of diffusion [59l p. 439-444]. Using 
a traveling wave ansatz u(x,t) = u(x — ct) =: u(£) 
for © yields the ODE 



du 



0. 



(7) 



Analyzing j7]) in the phase space variables (u, u') —: 
(u, v) shows that the point (u, v) = (1, 0) is a saddle 
and (it, v) = (0, 0) is a stable node or spiral. It is 
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straightforward 59, p. 441-442] to check that hetero- 
clinic orbits from (1,0) to (0,0) with u > 0, which 
correspond to non-negative traveling waves, can only 
exist in the stable node case for wave speeds c > 2. 
Since the orbit is directed from u = 1 to u = one 
also refers to this situation as the stable state u = 1 
invading the unstable state u = 0; see Figure QJa). 

Remark: Wave speeds c > correspond to waves 
traveling to the right. However, the FKPP equation 
(J6j> is invariant under the symmetry x — > — x so that 
a localized initial condition near x — triggers a pair 
of fronts, one traveling to the left one to the right. 

It is known that the traveling wave solutions to 
([6|) form the important solution set [30] Thm 1.4- 
1.5]. The minimal wave speed c* FKPP — 2 is the 
asymptotic speed of propagation for the FKPP non- 
linearity [3]- Since the wave speed is determined by 
the linearized problem 



du 
~dt 



a u 
dx 2 



(8) 



at the leading edge near (u, v) = (0, 0) as detailed 
in [81] p. 38-42] one also refers to the traveling wave 
with c* FK pp = 2 as a pulled front. For a more general 
nonlinearity f(u; n) pushed fronts can exist where 
the asymptotic wave speed c is larger than the linear 
spreading speed c* [81] p. 56]. The wave speed for 
both types is asymptotic and only achieved after a 
transient period. For pulled fronts the asymptotic 
expansion yields [HU p. 78] 



c(t) 



ki 
t 



_fc2_ 

t 3/2 



as t — > oo 



with explicitly computable positive constants fc^a > 
0. Hence, the wave speed is approached from below 
by a power law for pulled fronts. For pushed fronts 
the convergence to the asymptotic speed is exponen- 
tially fast [81] p. 74]. Another correction occurs when 
a cutoff for the reaction term is introduced [7] which 
leads to a logarithmic correction term. Furthermore, 
if an initial condition does not decay fast enough 
as \x\ — > oo then faster speeds than c* occur [81] 
p. 46]. In particular, for an initial condition decay- 
ing like 0{e~ a \ x \) for a > the speed increases as 
c(a) = 0(1/ a) as a -> [70]. 

The results for invasion fronts of the FKPP equa- 
tion already indicate that the variety of scaling be- 
haviors could be ideal to determine early- warnings. 
In fact, wave spreading in the stochastic case is even 
more intricate [48 ] I49 ] . 



4 Background - Stochastic PDEs 

As a stochastic generalization of (O the intuitive 
idea is to consider the equation 



du r^d 2 u , . , 

~di = dx 2 + ^^ U ' + ff ( u M x ' 



(9) 



where rj(x, t) formally represents the 'noise'. Here we 
consider two choices for the term rj(x,t). The sim- 
plest is to consider a real- valued (ID) Brownian mo- 
tion B(t) PS Chapter 8] with mean E[B(t)} = and 
covariance E[B(t)B(s)} = mm(t, s) for < s < t. 
Then white noise can be defined via r](x, t) = r](t) = 
B where the derivative is with respect to time and 
interpreted in the generalized sense [3] p. 52-53]. The 
covariance is E[r](t)r)(s)} — S(t — s) and one may then 
write ([9]) in two equivalent forms 



dt 

du = 



Dpi 



D 



i')x- 



WH - 



g(u)B(t), 
dt + g(u)dB. 



(10) 



The existence and regularity theory of pH|) is well 
understood [23] . If the noise should depend on space 
and time the theory is substantially more involved. 

One possibility is to consider a Hilbert space U 
(e.g. L 2 (M)) and a symmetric non-negative linear 
operator Q acting on U and define a U- valued Q- 
Wiener process W(t). If Tr(Q) < +oo there exists a 
complete orthonormal system {fk}kLi such that 



Qfk = Afe/fc, 



for k G N 



where {\k}kLi is a nonnegative bounded sequence. 
Then one may use the convergent sequence 



W(t) 



E 



\hBk(t)fk 



(11) 



with independent one-dimensional Brownian motions 
Bk(t) as a definition [66] p. 86-89]. As expected one 
has E[W(t)} = and E[W(t)W(s)] = mm(t, s)Q so 
that Q can be viewed as the covariance operator. In 
this case one may formally write © as looking for 
u = u(-,t) in the form 



= £>§-£ + nf(u) + g(u)W 



dt 

du = 



D 



d 2 u 
dx 2 



Hf(u) dt + g{u)dW. 



(12) 



which can be interpreted in a precise integral form 
[66l Section 5.1]. A well-developed existence theory 
for jnj} is available [H Thm 7.4, Thm 7.6]. 

One is tempted to take Q = Id to mirror the 
finite-dimensional case to obtain a 'white noise' pro- 
cess. However, Q = Id is not of trace-class (since 
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Tr(Q) = +00) and the series (fTTI) does not converge. 
However, one may construct a cylindrical Wiener 
process for Q — Id [66] p. 96-99] and characterize 
space-time white noise as W = rj(x,t) which has 
covariance lbL[r](x, t)ij(y, s)] = S(x — y)5(t — s). The 
existence and regularity theory for space-time white 
noise is slightly more involved and already leads to 
problems if x € M 2 [TO] p. 54]. Since we exclusively 
restrict toi£l these problems do not arise here 
and the existence theory works [2§J Section 6.1]. 

Remark: Instead of viewing the equation on func- 
tion spaces one may also consider an approach |84) 
where the solution u = u(x, i) is a real-valued ran- 
dom field which is a basically equivalent [39] ap- 
proach. In this case, one has E[W (x , t)W (y , s)] = 
min(i, s) min(:r, y) and that space-time white noise 
is 



ijxdl 



W(x,t) = v(x,t). 



5 Background - Numerical SPDEs 

First, we briefly review basic methods to solve the 
SPDE ([12"]) numerically for space-time white noise. 
The case (TIT)]) will follow as a special case. A natural 
first step is to start with a spatial discretization |37j . 
Consider a finite interval [x\, xn] C K for some N > 
1 and augment (|12[) with zero, reflective or periodic 
boundary conditions. Define (Ax) := (xn — xi)/N 
and consider the numerical solution Uj(t) ~ u(x\ + 
(j — l)(Ax),t). Then the space-discrete version of 
(fT2|) is a system of stochastic ordinary differential 
equations (SODEs) 



D N 



dt 




(13) 



for j = 1,2, ... ,N where {Bj}^ =1 are independent 
one-dimensional Brownian motions and the N x N 
matrix L depends on the boundary conditions. For 
reflection conditions [24] it follows that Ljj = — 2 if 
j G {2, 3, . . . , N - 1}, Ln = -1 = L NN , Uj = 1 
1 and Li 



if \i — j\ — 1 and L^ = otherwise. For periodic 



conditions u4\ one uses L, 



-2 for all j, L\ 



1 = Ljvij Lij — 1 if = 1 and Li 



N — 

otherwise. 



For zero boundary conditions 



the values it? 







are fixed, the first and last equation in (fT3 



obeys Ljj = -2 for j € {2, 3, . . . , N - 1}, Ly = 1 if 
\i — j\ = 1 and Lij = otherwise. For the simpler 
case (|T0|) one has a single Brownian motion so that 
Bj = B for all j and the factor 1/v 'Ax in (|13l) is 
removed [24] . 

It remains to solve the SODE flT3l which can be 
more compactly written as 



dU = F(U)dt + G{U)dB 



(14) 



where we view F(U) = (Fi(U), . . .,F N (U)) T , B = 
(dBi, . . . , dB N ) T as (column) vectors and G(U) is 
a diagonal matrix with diagonal entries Gjj(U). As 
a numerical scheme we shall always use either use 
the Euler-Maruyama method [35] or the the Milstein 
method in its explicit [42] p. 345-351] or implicit [42j 
p. 399-404] form stated below. 

Remark: The Milstein method is usually good 
as an exploratory tool due to its robustness. It has 
strong order-one convergence [42J, Thm 10.3.5]. It 
is relatively straightforward to implement the Mil- 
stein method as no multiple stochastic integral eval- 
uations occur since x € M [42, Chapter 10-11] [551 
p. 2]. Furthermore, it has recently been shown that 
it nicely extends to multiplicative trace-class noise 
[38] . Hence, the Milstein method provides a quite 
remarkable compromise between theoretical error es- 
timates and practical implementation issues; see also 
[32] for a computational introduction with test codes 
for scalar problems. The Euler-Maruyama method is 
faster but not as robust so it complements Milstein- 
methods nicely if many sample paths have to be cal- 
culated for a well-understood parameter regime. 

To state both schemes consider t £ [0, T] and de- 
fine (At) := T/K for some fixed K € N. Denote 
the numerical solution by U k ~ U(k(At)) where k £ 
{0, 1,2,..., K} and let (AB k ) = B((k + l)(At)) ~ 
B(k(At)) denote a vector of J\f(0, At) normally dis- 
tributed independent increments used at the k-th 
time step. The explicit Euler-Maruyama method is 
given by 



fe+i 



= XJ k j + At F 3 (U k ) + G. ]:j (U k ) (AB k )j (15) 

for each component j E {1,2,..., N}. The explicit 
Milstein scheme for (fT4]) - with diagonal noise matrix 
G - is given by \\j p.348,(3.12)] 



fc+i 



m + AtFj(U k ) + Gjj(U k ) (AB k )j 



+ \G n (U k ) d -§}f(U k )[(AB k ) 2 -At} 



(16) 



are discarded and the (N — 2) x (N — 2) matrix L 



For non-diagonal noise satisfying a suitable commu- 
tativity condition the scheme is still quite simple [42] 
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p.348,(3.16)] while for more general cases one has to 
be careful [38] . The implicit version of the Milstcin 
scheme for our problem is [42J p. 400] 

U k+1 = U k + At F 3 (U k +i) + Gn (□) (AB% 

+ \G 33 (U) a -§ j f(a)[(AB k ) 2 -At] [U) 

where we have the choice to make the scheme fully- 
implicit with □ = U k+1 or semi-implicit with □ = 
U k . Since the deterministic drift term F causes the 
stability problems if D(At) > (Ax) 2 24, p.64,67] it 
makes sense to chose the semi-implicit version. The 
algebraic problem to solve for U k+1 in ([T7|) can be 
solved using standard techniques such as Newton's 
method. 

It should be noted that the convergence and er- 
ror estimate of the numerical scheme do not imme- 
diately yield error estimates for quantitative proper- 
ties or scaling laws of traveling waves for the FKPP 
equation. For example, it has been demonstrated [l8j 
p. 71] that for pulled fronts of a discretized determin- 
istic FKPP (D = 1 = n) the speed given to leading- 
order by 

c* = 2-2(At) + ^(Ax) 2 + ■■■ . (18) 

A similar effect is expected for the stochastic FKPP 
equation and properties such as the diffusion proper- 
ties of the wave speed. Therefore, we have to view nu- 
merical scaling laws as approximations which carry 
some information about the discretization step sizes 
Ax and At. To minimize this effect, the formula ([T8]) , 
the stability requirement At < (Ax) 2 and the goal to 
minimize computation time indicate that we should 
choose At only slightly smaller than (Ax) 2 . 

We are not interested in computing the exact 
wave speed, only its trend under parameter variation 
will be relevant here. A simple method to compute 
an upper bound c on the wave speed for the initial 
condition u(0, 0) = 1 and u(x, 0) = for a stochastic 
wave is to collect the set of points (j Ax, k At) such 
that u(j Ax, k At) is less and u((j — 1) Ax, k At) is 
bigger than a threshold (usually we pick the thresh- 
old as 0.05). For each point one computes the esti- 
mate c w Ax j At and obtains c as the maximum. 

6 Linear Multiplicative Noise 

The first stochastic version of the FKPP equation we 
consider was studied by Elworthy, Zhao and Gaines 
[l"9 p 4" ] and is given by 

du fj, 2 d 2 u 1 . . „ • 

m=To^ + ^ uil - u) + euB (19) 



with one-dimensional time-dependent white noise B, 
a small parameter < fi <C 1 and noise strength 
e > 0. The multiplicative noise can be motivated 
e.g. by the interaction of a population u with the 
environment |78j or by parameter noise in the deter- 
ministic part [681 p. 8] such as a fluctuating growth 
rate [S3]. A multiplicative noise term of the form 
g(u) — u has also been used in a model for plank- 
ton spreading [S5] eq. (4b)] . For further mathematical 
background on SPDEs of the form (TT9"]) we refer to 
Section [4] 

It is proven in [19) that there are three major 
regimes for (|19[) depending upon the noise strength 
parameter e. For some n = 0(1) as /x — > 0, the cases 
e ~ k/u 2 , e ~ k/ii and e ~ k are identified as the 
strong, mild and weak noise regimes respectively. 
Elworthy, Zhao and Gaines prove and numerically 
demonstrate that for weak noise the wave propaga- 
tion of the pulled front is basically unaffected while 
the wave fails to propagate in the strong noise regime 
[T9l Thms 8.1-8.3]. In the mild noise regime the wave 
speed is decreased as \J2 — k [24] p. 65]. It is impor- 
tant to note that the spontaneous collapse of newly- 
introduced alien populations, which can occur in a 
strong noise regime, has been considered from an 
applied perspective in [76] . 

Using the scaling law of Brownian motion and 
the transformation 

x/i 2 2 
x i-> —=, t\->tfi , e := fie, 
V 2 

as discussed in Section [3j in the SPDE ([T9l) yields 
the more familiar form of the FKPP equation 

du d 2 u , . ■ , , 

- = — + u(l-u) + euB. (20) 

This gives the quite natural view that e>l, e~l 
and e <C 1 are the strong, mild and weak noise 
regimes. Figure[T]shows typical solutions for the three 
regimes; for details on the numerical methods see 
Section [5] The initial condition is taken as the intro- 
duction of a species at a particular fixed location so 
that 

u(x,t = 0) = H if f ^ = ' (21) 
y otherwise. 

It is understood that the numerical initial condition 
is obtained by choosing a mesh having a mesh point 
x = with u(0,0) = 1. 

Now consider the situation of the 'ecological ob- 
server' who can only measure the invasion wave at 
one particular point in space. Based on the results by 
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Fig. 1 Simulation of (f20)l using the implicit Milstein scheme (fl7T) with parameters K = 100, T = 20, N = 10 3 on 
the interval [—50, 50] with Neumann boundary conditions and initial condition u(x, 0) = 1 if x = and u(x, 0) = 
otherwise, (a) e = 0.02, (b) e = 0.3 and (c) e = 1.2. 



Elwortliy, Zhao and Gaines [151121] on propagation 
failure of the wave with increasing noise strength it 
is intuitive that the local statistics recorded at a fixed 
point carry information about the traveling front. 

For convenience we pick the point as x = 0. Con- 
sider a single sample path u(0,t). Let T denote the 
final time and consider the two basic statistics 



T— 1 Jt 



T—to Jto 



J to u(0,t)dt, 
!l[u{Q,t)-ufdt 



1/2 



Figure [5] shows an average of u over 200 sample paths 
which we denote by U. From the results it is becom- 
ing clear, once one compares the U plot with the 
c plot, that a decreasing mean population size does 
provide the expected early-warning sign for a de- 
creased invasion front speed. We also simulated the 
same case shown in Figure [2] for space-time white 
noise W in (f20"|) . The results are qualitatively sim- 
ilar with a slight quantitative shift towards faster 
waves at comparable noise strength. 

In both cases a relevant new result is that also 
the local fluctuations captured by the variance show 
a quite interesting behavior. Consider the scenario 
where the actual carrying capacity for the popula- 
tion is unknown. In this case, the population level 
u is insufficient to determine how far we are from 
propagation failure of the wave. Naively, one may 
interpret small population fluctuations as an indica- 
tor for a fast propagating wave but Figure [2] shows 
that it could equally well be a very slow propagating 
wave for a low population level. Hence one has to 
increase or decrease the noise strength to probe to 
which part of Figure [2] the observations match. 
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Fig. 2 Dependence of the time average U and the wave 
speed c on the noise strength e averaged over 200 sample 
paths. The SPDE (I20|l has been numerically solved (using 
Euler-Maruyama fig]!) with K = 100, T = 20, N = 10 3 on 
the interval [—50, 50] with Neumann boundary conditions 
and initial condition u(x, 0) = 1 if x = and u(x, 0) = 
otherwise. The top part shows U (circles) which has been 
calculated as the mean of the time series u(0, t) recorded 
by an ecological observer at the origin for t 6 [10,20]. 
The dots indicate ±1 standard deviation S for the time 
series; the curves are associated interpolations forming a 
confidence neighborhood. The bottom part of the figure 
shows an (upper bound) estimate for the wave speed. 



7 Nonlinear Multiplicative Noise 

As pointed out at the beginning of Section HI the 
noise terms euB and euW could be interpreted as 
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parametric or environmental noise. Another possible 
source of noise is 'individual-based' or 'finite-system- 
size' which we shall focus on in this section. Miillcr 
and Tribe showed in [S7J that the SPDE 



du 1 d 2 u 



u) + V2u W 



(22) 



arises as a limit of a contact process on a lattice 
originally studied in [5J as a model for long-range 
offspring displacement; traveling wave solutions to 
(f2"2"| exist for suitable parameter values 80 . How- 
ever, Miiller and Tribe also studied the behavior of 
(|22| with a scaled noise term ^Ju W varying the pa- 
rameter [i and proved [SHI Thm 1] that there exists 
a critical value /i c , independent of u(x, 0), such that 



(u(x, t) survives) = if /J, < fi c 
{u(x, t) survives) = 1 if fj, > fi c 



(23) 



Hence propagation failure of waves can occur like in 
the situation with noise term uW. Using the map- 
ping 1 i— > t//j,, x i — y x/y/Wfl, a standard scaling law re- 
sult [53 Lem 2.1.1] and the definition e := \/2(6/^) 1/4 
transform (12211 to 



du 
~di 



dx 2 



+ /m(1 - u) + e-y 



W. 



(24) 



Figure [3] shows the dependence of the population 
level, its fluctuations and the wave speed on the pa- 
rameter e. The results are very similar to Figure [2] 
with propagation failure for higher noise level as ex- 
pected from (|2"5|) . In particular, the conclusions from 
Section [HI about inferring wave propagation proper- 
ties from local data still apply. One may conjecture 
that the conclusions might apply to even more gen- 
eral versions of the FKPP equation with a noise term 
eg(u)W (or eg(u)B) as long as g(0) = 0. 



However, the SPDEs pj and (JUj) have noise 
terms that increase monotonically with the popu- 
lation level. This may not be realistic in all situa- 
tions as one expects the noise to change as u ap- 
proaches the carrying capacity. This is one motiva- 
tion to study the SPDE 



du d 2 u 



t) + ey/u(l - u) W. 



(25) 



This model was studied by several groups. In [55) 
it was proved for sufficiently small noise that com- 
pactly supported initial data remain within a time- 
dependent interval and that a well-defined front as 
well as an asymptotic wave speed exist; interestingly, 
the shape and asymptotic form of waves has been 
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Fig. 3 Dependence of the time average U and the wave 
speed on the noise strength e for J24j). Parameter values 
are as for Figure [2] 



of interest by an independent group for a discrete 
stochastic model [50] . 

Detailed numerical studies of the wave speed for 
(f2"5|) have been carried out [S21I5] focusing on the 
small noise regime and the fluctuation properties of 
the front. Due to the special structure of the FKPP 
equation one may also exploit a duality argument 
of (|25|) to a particle system [T5 1 I75 ] . Doering, Miiller 
and Smereka conjecture [15[ eq (55)] from the dual- 
ity relation that the asymptotic wave speed in the 
strong noise regime is given by c ~ 2/e 2 as e — > oo. 
Although it is unclear whether this conjecture is cor- 
rect it is evident from numerical simulations |15[ Fig 
2] that the wave speed decreases upon increasing the 
noise. 

Figure [4] shows the mean, standard deviation and 
wave speed calculated for (I25| . There are some mi- 
nor differences between this case and g(u) = u and 
g(u) — y/u shown in Figures [5] and [3] There is a 
larger plateau for small noise and it takes larger 
noise strengths to reach the vicinity of propagation 
failure. However, the main warning-signs from local 
data still remain as it is still possible to conclude 
from large population levels and low fluctuations a 
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Fig. 4 Dependence of the time average U and the wave 
speed on the noise strength e for (12511 . Parameter values 
are as for Figure[2]except for the slightly smaller time-step 
size N = 3 ■ 10 3 . 



fast wave while increasing fluctuations lead to slower 
speeds and low population levels with smaller fluc- 
tuations indicate closeness to propagation failure. 

Another different form of the noise term given by 
g(u) = u(l — u) was considered in [68] eq (39)-(41)] 
but we shall not consider it here as the results are 
similar. 

In summary, one should always measure the close- 
ness to carrying capacity and the size of the fluctu- 
ations (standard deviation, variance). If system pa- 
rameters change slowly one may determine from Fig- 
ures [20 whether the distance to propagation failure 
has increased or decreased. 



8 Transitions for the Allee Effect 

Although propagation failure is extremely interest- 
ing from the viewpoint of critical transitions, it is 
certainly not the only invasion wave phenomenon 
where local early-warning signs are desirable. As al- 
ready discussed in Section|3]there can also be pushed 
fronts if the nonlinearity f(u) is chosen differently. 
A reasonable prototypical model to study is the fol- 



lowing SPDE 

du d 2 u 



u){u — fi) + eg(u) W 



(26) 



which has been considered in |68j . The nonlinearity 
/(it) = u(l — u)(u + ji) may obviously arise due to 
an Allee effect in the context of ecology but it is also 
commonly used in other areas of mathematical biol- 
ogy, e.g. in neuroscience ([26]) would be referred to as 
Nagumo's equation (60]. We briefly recall some re- 
sults about the deterministic PDE (e = 0) described 
in [65]. For /i e [—1/2, 1/2] there exists a closed- form 
wave 

1 

u(x. t) — u(x — ct) = u(£) 



1 + e^ 



(«-£-) 



for an arbitrary phase £ > and propagation speed 



(27) 



There are three interesting special points. For fi — 
1/2 the front speed vanishes and the front reverses 
direction if /i > 1/2. The regime for /i e (0,1/2) 
is called bistable and changes to a pushed front at 
fi = 0; see [65] Sec 1] and references therein for 
details. The pushed front regime applies for \i e 
(—1/2,0) and at /i = —1/2 there is a pushed-to- 
pulled front transition [68] Sec V]. The wave speed 
for the pulled front is c* = 2^=71 [H Sec 3.3] gH 
p. 38-42]. Therefore, it is interesting to try to find 
early-warning signs for approaching the three spe- 
cial points /i = —1/2,0, 1/2 cases. The front rever- 
sal case fj, = 1/2 is clearly important as a direction 
change for an invasion front could be regarded as a 
critical transition but the other two cases could be 
of interest as well. 

For the SPDE ([26]) we shall choose the simple 
multiplicative noise g{u) — u. We consider the pushed- 
to-pulled transition first and try to apply our ap- 
proach from Sections EE Figure [5] shows the analog 
of the top parts of Figures [2][4] We observe that it 
is impossible to detect a trend or infer the speed of 
the wave. Hence, for the pushed-to-pulled transition 
with space-time white noise the classical variance- 
based early-warning signs cannot be applied for local 
data perturbed by a fixed noise level and observed 
at the center of the wave. Therefore, one should also 
think of new early-warning sign techniques in the 
context of wave propagation. 
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Fig. 5 Dependence of the time average U on the noise 
strength e averaged over 200 sample paths. The SPDE l|26l) 
for g(u) = u has been numerically solved (using Euler- 
Maruyama Jl5j) with K = 100, T = 15, N = 10 3 on the 
interval [—50, 50] with Neumann boundary conditions and 
initial condition u(x, 0) = 1 if x 6 [—1, 1] and u{x, 0) = 
otherwise. U (circles) has been calculated as the mean of 
the time series u(0, t) recorded by an ecological observer at 
the origin for t £ [7.5, 15]. The dots indicate ±1 standard 
deviation £ for the time series; the curves are associated 
interpolations forming a confidence neighborhood. 



Note carefully that we always used in our compu- 
tations in Figures [2][5] the regime for u(0, t) when the 
wave is already fully formed with t € [To, T] for some 
T 3> 1. However, the transient regime starting from 
the localized initial condition may also contain im- 
portant information. Figure [6] shows three numerical 
simulations for /1 = —0.3,0.2,0.4. 

The computation suggests that the initial tran- 
sient spreading of the wave u(0,t) for t £ [0,T ] is 
interesting. A simple measure to consider is 



/,„ := min {it(0, t) 

*e[o,T ] 



for a given u(x, 0)}. 



(28) 



Clearly, the result depends upon the choice of T 
and the initial condition u(x,0). However, if both 
are fixed then we may compare the results. Figure 
[7ta) shows the results for a parametric study of /j, £ 
[—0.75,0.5]. The insets (b)-(c) show a finer mesh 
resolution near the pushed-to-pulled transition at 
/1 = —1/2 and near propagation failure which is 
slightly shifted from the theoretical value at /1 = 1/2 
as the small hnite-width initial condition and the 
noise both seem to contribute to reach the absorb- 
ing state u = for parameter values fi smaller than 
1/2. In fact, due to these effects, the transition is 
more drastic than the formula (|27|) predicts. 

From Figure [7]Jb) it is apparent that the pushed- 
to-pulled transition is probably unpredictable from 
local data collected at x = 0. Since the wave speed 
transition is continuous one should probably not clas- 
sify the pushed-to-pulled transition as a 'critical tran- 
sition'. Therefore, it is not crucial to predict it but 
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Fig. 7 (a) Dependence of the minimum u m defined in 
1281 on with To = 10, averaged over 200 sample paths. 
The SPDE (I26|l for g(u) = u has been numerically solved 
(using Euler-Maruyama {HJ) with K = 100, T = 20, 
e = 0.05, N = 10 3 on the interval [-50, 50] with Neumann 
boundary conditions and initial condition u(x, 0) = 1 if 
x £ [—1, 1] and u(x, 0) = otherwise. The circles indicate 
u m and the dots ±1 standard deviation £ calculated from 
the sample paths, (b) Zoom near the theoretical pushed- 
to-pulled transition at fi = 1/2. (c) Zoom near propaga- 
tion failure transition. 



the result shows the limitation of the ecological ob- 
server at x = 0. The same conclusion applies to the 
change from the pushed to the bistable regime at 
/i = 0. For the propagation failure scenario Figure 
[Tie) shows a scaling law for the decrease of u m and a 
slightly increasing variance may help us to anticipate 
the upcoming critical transition. This should not be 
surprising since we already considered similar prop- 
agation failure cases in Sections EE The difference 
is that we used a completely different indicator in 
Figure 0[c). 

As for the classical FKPP equation one should re- 
mark for the Allee effect situation that different noise 
terms certainly do make sense, e.g. g(u) = u(l — u) 
considered in [2j eq.(2)]. Based on the observations 
for varying the noise terms for the classical FKPP 
equation, and obtaining similar results for several 
choices, we shall not consider these generalizations 
for 



9 Noncompact Initial Invasions 

Based on the results in Section [S] we have observed 
that the initial transient regime, starting from a lo- 
calized invasion wave, can be useful. It remains to 
consider the case when the initial condition is not 
localized. In particular, we consider the FKPP equa- 
tion 



du d 2 u 

m^dx^ + u{1 



euW. 



(29) 
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Fig. 6 Simulation of l|26l) for g(u) = u using the implicit Milstein scheme (11711 with parameters K = 100, T = 20, 
iV = 2 • 10 3 on the interval [—50, 50] with Neumann boundary conditions and initial condition u(x,0) = 1 if x € [—1, 1] 
and u(x, 0) = otherwise, (a) fi = —0.3, (b) fi = 0.2 and (c) /i = 0.4. 
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Fig. 8 Dependence of the minimum u m defined in !28l on a 
for (|29II - (I30|| ; average over 200 sample paths. The param- 
eters for the numerical simulation(using Euler-Maruyama 
H15jl ) on the interval [—50,50] with Neumann boundary 
conditions are K = 150, T = 10, e = 0.05 and N = 500. 
(a) The circles indicate u m and the dots ±1 standard de- 
viation calculated from the sample paths, (b) Wave speed 
c (circles) and associated ±1 standard deviation (dots). 



with initial condition 

u(x,0) = e~ alxl , farar>0. 



(30) 



Recall from Section [3] that the wave speed scales as 
c{a) = O(l/o) for a -> 0. 

Figure EJa) shows the dependence of the initial 
transient observed at x = on the complete ini- 
tial data. Therefore, small minimum values indicate 
comparatively slower waves and no response to the 
initial condition (u rn m 1) signals a very fast wave. 



Figure|SJb) shows an upper bound to the wave speed 
and raises the interesting question whether we should, 
or should not, view a blow-up point for the wave 
speed as a critical transition. 

There are two main conclusions from the results 
for (J29j)-((30j) and from Section U Firstly, one defi- 
nitely should try to measure an invasion wave im- 
mediately once the first occurrence of a new popula- 
tion in a new environment has been observed. Sec- 
ondly, knowing the basic structure of the initial con- 
dition can be crucial for prediction, e.g. in Figure [7] 
<C u m < 1 still indicates a well-defined asymptotic 
wave speed in the pushed regime while for Figure 
[8] the condition <C u m < 1 indicates closeness to 
a wave speed blow-up point. Hence, it is crucial to 
know, on a qualitative level, whether the initial in- 
vasion is really localized or whether it really consists 
of a full front. 



10 Outlook 

Since early-warning signs and stochastic scaling laws 
for noisy traveling waves are still a relatively new di- 
rection, we have only been able to cover a few aspects 
here. Many open problems arose which we summa- 
rize here. 

The restrictions to one spatial dimension x € K 
and one population component u have to removed in 
the future. There are many interesting cases e.g. multi- 
component systems such as reaction-diffusion mod- 
els with predation [51] , FKPP-type plankton dynam- 
ics [6] or Nagumo (Allee effect )-type equations [2 7) . 
Multiple spatial dimensions can lead to more com- 
plicated bifurcation structures [40] . One may also re- 
move all restrictions which can generate interesting 
life-death transitions for multi-component, 2D and 
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3D systems [SS] . Another highly relevant generaliza- 
tion are heterogeneous [69] and random [85] envi- 
ronments. Furthermore, the structure of the FKPP 
equation may be too restrictive which suggests to 
add transport/advection terms and active bound- 
aries in which case discontinuous wave speed transi- 
tions have been reported pj] . Also the assumption of 
time-white or space-time-white noise is too restric- 
tive and one should extend the view to spatially- 
colored noise [25] and trace-class covariance oper- 
ators. Another issue that looked interesting is the 
relevance of fluctuations ('front diffusion') [TIIBT] for 
early-warning signs. 

In all cases, our main driving question in this 
paper seems to be open: How much information do 
local statistics of an SPDE, collected at one (or mul- 
tiple) locations, carry about the speed and bifurca- 
tions of traveling waves? It is seems plausible to ob- 
tain basic answers to these questions using numer- 
ical simulations. To develop a mathematical theory 
for quantitative scaling laws of SPDEs and their ap- 
plication to critical transitions is expected to be a 
challenging problem for a long time. 
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